Draft version August 21, 2012 

Preprint typeset using L^T^X style emulatcapj v. 5/2/11 



CAN PLANETARY INSTABILITY EXPLAIN THE KEPLER DICHOTOMY? 



(N 

o 

(N 

< 
00 



PL, 

Of 

6 

o3 



(N 
> 

oo 
On 1 
00 

O 

(N 



13 



Anders Johansen, Melvyn B. Davies, Ross P. Church, and Viktor Holmelin 

Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 221 00 Lund, Sweden 

Draft version August 21, 2012 

ABSTRACT 

The planet candidates discovered by the Kepler mission provide a rich sample to constrain the ar- 
chitectures and relative inclinations of planetary systems within approximately 0.5 AU of their host 
stars. We use the triple-transit systems from the Kepler 16-months data as templates for physi- 
cal triple-planet systems and perform synthetic transit observations, varying the internal inclination 
variation of the orbits. We find that all the Kepler triple-transit and double-transit systems can be 
produced from the triple-planet templates, given a low mutual inclination of around five degrees. Our 
analysis shows that the Kepler data contains a population of planets larger than four Earth radii in 
single-transit systems that can not arise from the triple-planet templates. We explore the hypothesis 
that high-mass counterparts of the triple-transit systems underwent dynamical instability to produce 
a population of massive double-planet systems of moderately high mutual inclination. We perform 
A^-body simulations of mass-boosted triple-planet systems and observe how the systems heat up and 
lose planets by planet-planet collisions, and less frequently by ejections or collisions with the star, 
yielding transits in agreement with the large planets in the Kepler single-transit systems. The result- 
ing population of massive double-planet systems can nevertheless not explain the additional excess 
of low-mass planets among the observed single-transit systems and the lack of gas-giant planets in 
double-transit and triple-transit systems. Planetary instability of systems of triple gas-giant planets 
can be behind part of the dichotomy between systems hosting one or more small planets and those 
hosting a single giant planet. The main part of the dichotomy, however, is more likely to have arisen 
already during planet formation when the formation, migration or scattering of a massive planet, 
triggered above a threshold metallicity, suppressed the formation of other planets in sub-AU orbits. 
Subject headings: planets and satellites: formation - planets and satellites: dynamical evolution and 
stability - protoplanetary disks 



1. INTRODUCTION 



Radial velocity surveys searching for exoplanets have 
been operating for more than 15 yea r s sinc e the first 
exoplanet discovery (jMavor fc Quelozl [1995!) . One of 
the most important discoveries of these surveys is that 
many stars host multiple planets within a few AU 



iButler et aHl!99ft IWright et alj|2009t lLovis eta! 
However, most of the detected planetary companions are 
massive since such planets produce higher radial velocity 
signals. Recent years have seen the discovery of a new 
class of low-mass exoplanets - super-Earths with masses 
less than ten times the Earth's - from radial v elocity sur- 
veys (e.g. lRivera et al.l2005HUdrv et al .12001 . Although 
the number of super-Earths dis covered by radial v elocity 
surveys is growing rapidly (e.g. lMavor et al.ll2011| ). their 
number is not yet high enough to extract statistical infor- 
mation about the abundance and general architectures of 
planetary systems including low-mass planets. 

The Kepler mission (which we will often refer to sim- 
ply as Kepler) provides the first statistically-significant 
survey of gas giants, Neptune-size planets . super-Earths 
and t errestrial planets orbiting other stars dBorucki et al.l 
I2010D . Monitoring 156,000 stars for periodic dips in the 
light curves, Kepler is sensitive to planets as small as one 
Earth radius (|Fressin et a,lJl20T2h . The relative numbers 
of systems with single and multiple transiting planets is 
a sensitive function of the intrinsic multiplicity of plane- 
tary systems and the relative inclinations between planet 
orbits and thus hold important information about plan- 



etary system architectures. 

The target stars of Kepler are typically too faint for 
follow-up observations of all but the largest planets, so 
detections by Kepler are generally referred to as planet 
candidates, unless their masses ca n be determined b 
radial velocity measurements fe. g. IBatalha et all 1201 
or transit timing variations (e.g. iLissauer et al.l 12011 ' 
However, due to the high photometric precision, fewer 
than 10% of the planetary candidates are b elieved to be 
false candidates (|Morton fc Johnson! 120 111) . The strong 
clustering of planet candidates into multiple systems is 
further evidence for the physical nature of the detec- 
tions, since false positives caused by binary stars would 
be distributed even ly among the Kepler target stars 
(|Lissauer et al.ll2012[) . 

The first four mont hs of Kepler data re vealed 1235 
planetary candidates (|Borucki et al.l l2011al ) while that 
number had grown to 2321 after sixteen months 
(jBatalha et al.N2012t l - 253 of Earth size (R < 1.25-R®), 
712 of super-Earth size (1.25-R® < R < 2R m ), 1078 
of Neptune size (2i?® < R < 6i?®), 207 of Jupiter 
size (6i?® < R < 15i?@) and 71 with sizes larger than 
Jupiter (15-R® < R < 22i?@). A large fraction of those 
planets are in mult iple systems (jLissauer eTaTll201ibl : 
IBatalha et al.ll20T2T ) - 245 host stars show double tran- 
sits, 84 show triple transits, 27 show quadruple transits, 
8 show quintup l e tran sits and 1 shows sextuple transits. 
lOfir fc Dreizlerl (|2012l ) used an alternative data reduc- 
tion algorithm to analyse the Kepler data and found sev- 
eral new transit signals, upgrading one of the planet host 
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stars to a new candidate sextuple system. 

The Kepler data has already led to new knowledge 
about the characteristics of the innermost parts of plan- 
etary systems. Individual systems such as the sextu- 
ple Kepler- 11, with masses deter mined by transit timing 
variations (jLissauer et al.ll2011aD . show that small plan- 
ets can have a variety of densities and hence represent 
distinct classes of planets. Statistical analysis of Kepler 
data gi ves cons t raints on formation and migration pro- 
cesses. lYoudinl (|2011f ) showed that there is a difference 
between the size distributions of shorter and longer pe- 
riod planets. The apparent lack of planets of around 
three Earth radii in size, relative to both smaller and 
larger planets, in orbits shorter than 7 days may be evi- 
dence of either sublimation of volatiles or inefficient gas 
accretion in Neptune-mass planets migrating to close or- 
bits. 

Transit observations are mainly biased by the orbital 
alignment with the line of sight (affecting all planets) and 
stellar and instrumental noise (affecting mainly smaller 
planets). The orbital alignment bias is relatively inde- 
pendent of planet size, since stellar radius and planetary 
semi-major axis dominate the transit probability. The 
transit bias is thus less complicated than the bias in ra- 
dial velocity detections. Together with the sensitivity 
down to Earth-sized planets, this makes the Kepler data 
well suited for studies of planetary system architectures. 

The goal of this paper is to find underlying planetary 
system architectures that explain the absolute and rela- 
tive number of single and multiple transits in the Kepler 
data. Our approach d iffers from that of other authors. 
iLissauer et al.l (|2011b|) used the Kepler data to param- 
eterise the population of planets in terms of their sizes 
and semi-major axes. They found that the majority of 
the observed double-transit systems and a substantial 
fraction of the observed single-transit systems could arise 
from relatively flat systems of higher multiplicity, while 
the remaining single -transit systems come from a second 
distinct population. iTremaine fc Donel ()2011l ) looked at 
the problem from a statistical approach. They concluded 
that the inversion from observed transits to underlying 
populations is in principle degenerate in the mutual incli- 
nation parameter, in that solutions consisting of a com- 
bination of single-planet systems and very densely popu- 
lated systems (up to 40 planets within 0.5 AU) of nearly 
isotropic orbits can yield the observed transit rates. It 
is nevertheless not clear that these nearly isotropic so- 
lutions arc physical, in that most planets need to be in 
very densely packed syste ms, with few or no plan etary 
systems of low multipli city. IWcissbein et all (120121) used 
an approach similar to ILissauer et al.l (|2011b| ) and find 
that their assumption that planetary occurrence is an 
independent statistical process is not supported by the 
data. In this paper we use an alternative approach to 
constraining underlying planetary system architectures. 
We use the triple-transit systems observed by the Ke- 
pler mission as templates for physical planetary systems 
and perform synthetic transit observations to derive the 
relative number of single, double and triple transits as a 
function of the mutual inclination between the planet or- 
bits. We thus work close to the original data and do not 
have to go through steps of interpretation and analysis 
of the statistical properties of planet sizes and orbits to 
fit these with distribution functions. 



Our paper is organized as follows. In <J2] wc explain 
how we select the underlying planet population by using 
the triple-transit systems observed by the Kepler mission 
as templates. In JJ] we perform synthetic transit obser- 
vations of these systems, weighting each template system 
by the inverse of the triple-transit probability. We com- 
pare the synthetic transits with the Kepler data in ^land 
find that there is a good match in radii and semi-major 
axes between the synthetic double and triple transits and 
the Kepler data. However, the synthetic single transits 
do not match those observed - the synthetic catalogue 
contains too few single transits (by about a factor of 
three) and it completely lacks planets having large radii 
(between the sizes of Neptune and Jupiter) which have 
been found by the Kepler mission. This suggests that in 
addition to the population of three-planet systems (seen 
in double and triple transits), there is a second, distinct 
population of planetary systems containing larger plan- 
ets. We discuss the properties of this second population 
in ij5] We explore the hypothesis that planetary sys- 
tems form with a range of masses and that the most 
massive are inherently unstable to planet-planet interac- 
tions. We perform iV-body simulations of mass-boosted 
versions of the Kepler triple-planet templates in SJ5] The 
typical result is that two planets collide and merge and 
leave a moderately inclined double-planet system. While 
this second population in principle can produce the large 
planets seen in single-transit systems in the Kepler data, 
this does not explain an additional excess of small planets 
(smaller than four Earth radii) among the single-transit 
systems, nor the lack of gas giants in double-transit and 
triple-transit systems. In fJ7] we relax the assumption 
that all planetary systems are triple, constructing single- 
planet and double-planet systems by removing one or 
two planets from the triple-planet templates, and show 
that this does not change the necessity for a dichotomy 
of planetary systems in nature to explain the data. Wc 
conclude in ££SJ proposing that the main part of the di- 
chotomy between systems showing single transits and 
systems showing double transits and triple transits arose 
already during planet formation, when the migration or 
formation of a large planet suppressed the formation of 
additional planets in sub-AU orbits. 

2. SELECTION OF UNDERLYING PLANET POPULATION 

We will in this paper make the assumption that an i- 
multiple system observed by the Kepler mission can be 
used as a template for a physical i-multiple system. A 
system with i transiting planets observed by the Kepler 
mission may in reality have more planets than the ob- 
served i, especially further from the star where the tran- 
sit probability is low, but a reduction in planet number 
produces a system that is equally physical, as such re- 
duction very rarely leads to dynamical instabilities. Wc 
sometimes use the notation it-jp for observing i transit- 
ing planets in an intrinsically j-multiple system. 

To get a uniform and unbiased sample of planetary 
systems, we first apply a number of selection criteria to 
the planets in the 1 6-months Kepler data: 

1. The signal-to-noise ratio must be larger than 16. 

2. The planet period must be shorter than 240 days. 

3. The planet radius must be smaller than 22.4 Earth 
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radii. 

The first criterion ensures that the sample is close to 
being complete and that no size bins are dominated 
by fortuitous detections, an important issue when com- 
paring observations of synt hetic planetary system s to 
the Kepler data. We follow iLissauer et al.l (|2011b[ ) and 
choose a limiting signal-to-noise ratio of 16. The second 
criterion includes only planets seen to transit at least 
twice in the 16-month data. The third criterion ensures 
that only planetary-mass objects are selected, applying 
a maximum radius o f 22.4 Earth radii, as proposed by 
IBorucki et all (|2011bf ). 

Selection criteria are applied on a planet-by-planet ba- 
sis, so dense systems may become sparser after the selec- 
tion criteria. The triple system KOT284 has two planets 
in very similar orbits. We exclude this system because 
the two planets likely orbit the two differ e nt stellar com- 
pone nts of a binary (jHowell et al.1 120111 : ILissauer et al.l 
l 2012t ). Another peculiar system is KOI-191, an originally 
quadruple system with a Jupiter-sized and an Earth- 
sized planet very close to each other. However, the small 
planet has a low signal-to-noise ratio and is excluded by 
criterion 1 abov e, reducing KOI - 191 t o a regular triple- 
transit system. iFabrycky et al.1 ([20121 ) reported a simi- 
larly curious quadruple system, KOT2248, with two very 
close planets, but the detections of two of the four plan- 
ets have a low signal-to-noise ratio, and hence KOT2248 
reduces to a regular double-transit system after applying 
our selection criteria. 

Fig. [T] shows the triple-transit, double-transit, and 
single-transit systems present in the Kepler data after ap- 
plying selection criteria. We do not consider higher-order 
systems since these are too few to serve as a statistically 
significant template base for the synthetic systems. The 
number of triple systems is reduced from originally 84 to 
62, the number of double systems is reduced from 245 to 
187, while the number of single systems is reduced from 
1425 to 1183. 

The orbital periods and radii of the selected planets are 
shown in Fig. [2] While planets in systems with two or 
three transits are statistically similar, planets in systems 
with only a single transit are clearly on average signif- 
icantly larger than t he planets in higher- order systems 
([Latham et al.ll20lH ILissauer et al.ll2011bf ). 

3. CALCULATION OF TRANSIT RATES 

The goal of the paper is to perform simulated transit 
observations of populations of synthetic planetary sys- 
tems to find the intrinsic population of planetary sys- 
tems that matches the Kepler data. We use the triple 
planetary systems observed by Kepler, after application 
of selection criteria to make the sample uniform (see Sj2]), 
as templates for physical systems and assume that stars 
have either zero planets or a planetary system with the 
radii and semi-major axes of one of the triple-planet sys- 
tems observed by Kepler. Although crude, the assump- 
tion of planetary systems having either zero or three 
planets turns out to give a very good fit to the ob- 
served double-transit systems. Our approach is also very 
straight-forward to apply to the data, in that we do not 
need to go through complex steps of extracting statistical 
information about planetary sizes and semi-major axes 
from the data. We will later relax the assumption that 



planetary systems are always triple and show that this 
does not change the conclusions (sec $7]). 

All the template planetary "systems" are limited to 
a maximum semi-major axis of approximately 0.5 AU. 
While this is much smaller than the expected total extent 
of planetary systems in general, the Kepler data shows a 
surprising population density and variety in these rel- 
atively close orbits. We choose triple-transit systems 
as templates, since denser systems (quadruple and be- 
yond) arc observed in very low numbers and arc hence 
much more affected by small number statistics than the 
triple-transit systems. A s discussed in the introduction, 
iTremaine fc Donel ([201 If ) showed that very densely pop- 
ulated systems of nearly isotropic orbits can in principle 
explain the Kepler data, but we will in this paper assume 
that systems more dense than three planets are rare and 
hence do not contribute significantly to the triple transits 
in the Kepler data. This assumption is discussed further 
in SJSJ 

The orbits of all three planets are initially put in the 
same plane. We assume that the orbits have very low 
eccentricities between and 0.01 (a good approximation 
to the sta tistical properties of t he Kepler transit dura- 
tions, see iMoorhead et al.ll201ll) . For each planet in a 
triple system we then choose an inclination angle dis- 
tributed evenly between and /3. Here /3 is a measure 
of the opening angle of the system. We incline the or- 
bital plane of the planet by this angle relative to the 
initial, common orbital plane. Afterwards we choose a 
random angle between and 2tt and rotate around the 
axis perpendicular to the original plane to get a random 
longitude of ascending node and argument of periapsis. 

Fig. [3] shows how a chosen j3 = 5° yields a distribution 
in relative inclinations between planet pairs. We have 
taken pairs of planetary orbits and inclined and rotated 
the planes as described above. The distribution peaks 
at a mutual inclination around i ~ 3.5°. The maximum 
relative inclination between any two planets, 2/3, is ob- 
tained by first inclining both planetary orbital planes by 
j3 and then choosing opposite longitudes of the ascending 
node. 

We calculate the probabilities for an observer to see 
each of the planets transiting by considering 1,000 ran- 
dom realisations of each triple-planet system, for a given 
value of (3, from a large number of random directions. 
We choose directions uniformly distributed on the unit 
sphere and calculate the total number of single, double 
and triple transits as well identifying which planets are 
involved. 

Using it systems as templates for ip systems requires 
us to take into account that the it system was observed 
at all. Take for example an observed system with three 
planets close to the star and another observed system 
with three planets far away from the star. These two 
transiting systems can be used as templates for two phys- 
ical systems. However, if the wide system is ten times 
less likely to be observed than the tight system, then we 
must convert the two observed transiting systems into 
one tight synthetic planetary system and ten wide syn- 
thetic planetary systems. 

We need to go through a number of steps from the 
transit probability of individual planets to calculating 
the number of single, double and triple transits that are 
observed. We use a notation defined as following. 
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• Probability that planet i transits in system j: 

Pij = — ■ 

• Probability of exactly i transits in system j: 



(1) 



Pij=PiM- (2) 
Mean weighted probability of an i- transit: 



Pi = 



Weight of system j: 



Wi 



1 

Pv 



(3) 



(4) 



Here R+j is the radius of the host star in system j and 
is the semi- major axis of planet i in system j, with i = 1 
denoting the inner planet, i = 2 the middle, and £ = 3 
the outer planet. The probability ptj that planet i is 
seen to transit in system j is then given by equation (fT]) , 
ignoring the small contribution of the planetary radius 
to the transit probability. The probability that system j 
shows exactly i planets transiting is denoted Py (equa- 
tion [2]). This probability is a (complicated) function of 
the mutual inclination parameter (3. So far we have not 
considered the relative weight of the systems, but as ex- 
plained above, the fact that the Kepler mission observes 
a given triple-transit system automatically requires that 
systems with low triple-transit probabilities are more fre- 
quent among the target stars than systems with high 
triple-transit probabilities. Hence we define Pj as the 
mean weighted probability that any given host star shows 
exactly i planets transiting, given by a straight-forward 
weighting in equation ([3]). The weighting function is 
defined in equation ((4]). This choice ensures that each 
triple-transit system in the Kepler data will be observed 
exactly once in the synthetic observations. Effectively 
Wj = l/Psj is the number of randomly oriented copies 
of system j among the Kepler target stars. We thus im- 
plicitly assume that higher-order planetary systems do 
not contribute significantly to the observed triple-transit 
systems. 

The probabilities of observing each of the Kepler triple- 
transit systems, P$j using the above notation, are shown 
in Fig. 21 The highest triple-transit probability - ap- 
proximately 11.2% at p = 0° and 10.1% at P = 5°, is 
obtained for the KOI-1835 system, with three close-in 
planets orbiting a relatively large host star of radius 1.66 
Rq. Another noteworthy system is KOI-961 which has 
three small planets very close to the star, with periods 
0.45 days, 1.22 days and 1.87 days. However, the prob- 
ability of Kepler seeing three transits from this system 
is not as high as one might expect, sin ce the radius of 
the h ost star is low (0.17 solar radii, see iMuir head et all 
120121 ). The contrast between the lowest and highest 
triple-transit probability is around 10 for a mutual incli- 
nation parameter of p = 0° and 250 for p = 5°. Planets 
in wide orbits are much more likely to be inclined out 
of view as j3 is increased, compared to planets in tight 
orbits. 



4. FITTING THE TRANSITS OF THE KEPLER MISSION 

The transit detection algorithm described in $3] allows 
us to calculate the number of single, double and triple 
transits arising from the triple-planet templates and to 
compare the results with the Kepler data. We vary the 
mutual inclination parameter f3 between and 10 de- 
grees, in steps of 1 degree. Each synthetic triple-planet 
system is given a weight equal to the inverse triple-transit 
probability at the given inclination parameter. The re- 
sulting numbers of single-transit and double-transit sys- 
tems produced are shown in Fig. [SJ We normalize the 
number of single-transit and double-transit systems by 
the number of triple-transit systems. Matching those 
ratios to the Kepler data, the correct absolute number 
of transits can be obtained by assuming that a frac- 
tion / = (Ns/P^/N-,. of the Kepler target stars host 
triple-planet systems. It is clear from Fig. [5] that we 
can not simultaneously reproduce both the number of 
double-transit and single-transit systems in the Kepler 
data. Reproducing the correct fraction of single-transit 
systems overproduces the number of double-transit sys- 
tems by a factor two, while a faithful reproduction of 
double-transit systems will only produce 1/3 of the ob- 
served single-transit systems. An important property of 
fitting synthetic transits to observed transits is that no 
system produces a negative number of transits. Hence 
no population of planets is allowed to overproduce the 
number of single, double or triple transits found in the 
Kepler data. That puts a hard limit of j3 < 5° to avoid 
overproduction of double-transit system^. 

The radius and semi-major axis distribution of the syn- 
thetic transits are compared to the observed transits in 
Fig. |6l We focus on a mutual inclination parameter of 
P = 5° which approximately matches the ratio of double- 
transit to triple-transit systems in the Kepler data. The 
radii of planets in our observed triple-transit systems in 
Fig. [6] match by definition the observed distribution of 
planet radii in the template triple-planet systems. A 
similarly good correspondence is found in the double- 
transit systems. This way the triple-planet systems can 
both produce the correct relative number of double-to- 
triple transits and the at the same time reproduce the size 
distribution of the planets in double-transit and triple- 
transit systems. 

However, planets in single-transit systems seen in the 
Kepler data include a much larger fraction of planets 
having larger radii (between the radii of Neptune and 
Jupiter). The synthetic systems almost entirely lack 
planets with such large sizes. In addition, the number 
of single-transit systems produced from the triple-planet 
templates is about a factor of three lower than that ob- 
served. A higher number of single-transit systems could 
in principle be produced by choosing a higher mutual in- 
clination parameter P sa 8° in Fig. [51 however this would 
not produce a good match to the observed distribution 
of planetary radii, as evidenced in Fig. |6l 

Fig. [6] also shows the distribution of semi- major axis of 

1 This limit relies on our assumption that planetary systems 
with intrinsic multiplicity larger than three are rare. Consider- 
ing te mplates of higher multip licity would lead to a higher value 
of j3 (Trcmainc & Dong 2011). However, this would be in some 
conflict with the measured coplanarity of systems such as the sex- 
tuple Kepler- 11 where the mutual inclination is around one degree 
ULissauer et alj|20lTah . 
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the synthetic planets and the Kepler data. Again there is 
a perfect match, by definition, for triple-transit systems, 
whereas single-transit and double-transit systems have 
planets that are on the average a bit further from the 
star in the Kepler data. This discrepancy is nevertheless 
minor. 

5. PROPERTIES OF A SECOND POPULATION 

The inability for the template triple-planet systems 
to explain simultaneously the double-transit and single- 
transit systems in the Kepler data implies that there is 
a dichotomy in the underlying planetary systems. On 
the one hand there is a population of triple-planet sys- 
tems with small planets which faithfully reproduces all 
triple-transit and double-transit systems, but only 1/3 of 
the single-transit systems. An additional second popula- 
tion of planetary systems, distinct from the Kepler triple- 
transit systems, produces the remaining single-transit 
systems. The cumulative distribution and probability 
distribution functions of the planetary radii for observed 
and synthetic single-transit systems, in terms of the ab- 
solute number of planets, are shown in Fig. [7] We choose 
a mutual inclination parameter of /3 = 5° which matches 
the number of double-transit systems in the Kepler data. 
The synthetic data misses around 750 single-transit sys- 
tems. Of these approximately 500 are small (smaller than 
four Earth radii) while 250 are large (larger than four 
Earth radii). 

More insight into the second population of planets 
needed to produce the missing single-transit systems 
is gained by considering the number of single-transit, 
double-transit and triple-transit systems as we increase 
the mutual inclination parameter p. We maintain the 
triple-planet systems as templates for illustration, al- 
though we show in the next section that unstable triple- 
planet systems evolve to stability by having two or more 
planets collide. 

In Fig. [5] we show the transit number for a mutual in- 
clination parameter up to 40° . We maintain the normal- 
isation that the 62 triple-transit systems found in the 
Kepler 16-months data after applying selection criteria 
must be produced at the mutual inclination parameter 
P' = 5°. The transit numbers Ni are thus defined as 



(5) 



with iV 3 (/3') = 62 to match Kepler observations. Here 
the Pi are the mean i-transit probabilities per system, 
with each system weighted by its triple-transit probabil- 
ity at (3 = p' . The normalisation with N3 and P3 at 
P = P' implies that the population of planetary systems 
with mutual inclination angle f3 formed with /3 = /3', just 
like the flat population discussed above, and then later 
heated up to have a higher mutual inclination parameter 
of P > P'. The normalisation further implies that the 
equation for the transit number can be simplified as 



(6) 



with the short hand notation A3 = N 3 (P') and Wj = 
Wj(P'). Fig. [8] shows that the number of triple-transit 
systems falls rapidly with increasing p and is unimpor- 



tant (less than 10% of the value at P = 0°) at mutual in- 
clinations above P « 12°. The number of double-transit 
systems falls much more slowly. 

In the isotropic limit where planetary orbits arc inde- 
pendent we have 
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These numbers can be calculated from the output of the 
synthetic transit observations. The isotropic limits are 
shown in dashed lines in Fig. [H In this limit there are 
937 single-transit systems, 34 double-transit systems and 
0.5 triple-transit systems. At p = 40° there are neverthe- 
less still 67 double-transit systems, while the number of 
single-transit systems has risen to almost 900. In order 
to produce the 750 missing single-transit systems from 
a population of triple-planet systems of high mutual in- 
clination, a mutual inclination parameter of at least 20° 
would be needed. This high opening angle would as a 
side effect produce more than 100 double-transit systems 
additional to the ones produced from the low-mass triple 
systems, enough to push the number of synthetic double- 
transit systems significantly beyond the observed num- 
ber. These extra double-transit systems would further- 
more have the characteristic large planetary radii of the 
single-transit systems, in conflict with the excellent ra- 
dius match that we find between observed double-transit 
systems and synthetic double-transits from the Kepler 
triple-systems. Hence it is very difficult for a population 
of mutually inclined, massive planets in triple systems to 
explain the observed dichotomy between systems show- 
ing one transit and systems showing two or three transits. 

In the following we explore the hypothesis that the sec- 
ond population constituted the upper end of a continuous 
distribution of planetary birth masses, and that these 
systems of massive planets became unstable to planet- 
planet interaction and lost one planet by collision or ejec- 
tion, turning into double-planet systems with a moder- 
ately high inclination between the planetary orbits. A 
higher mutual inclination will favour observations of sin- 
gle transits over observations of double transits, and this 
way a second population of double-planet systems can 
contribute to the missing single-transit systems without 
polluting the double-transit systems which are already 
produced from low-mass triple-planet systems. 

For our hypothesis it is not important whether the 
massive planets formed in situ close to their host stars 
or whether they formed further out and migrated in due 
to gravitational torques from the gaseous protoplanetary 
disk. The nebula from which the solar system plan- 
ets formed did not contain enough mass close to the 
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star to form th e solid cores of multiple gas giants there 
(Hayashi 1981), although this property is a direct conse- 
quence of the observed lack of massive planets in sub-AU 
orbits in the solar system. Extrasolar protoplanetary 
disks are observed to contain a wide range of masses, 
up to 10 times higher than the minimum mass solar 
nebula, although such mas sive disks are relatively rare 
(| Andrews fe Williamsl [2005). Migration of massive plan- 
ets o ccurs on the visc ous time-scale of the protoplanetary 
disk (|Lin e t al. 199(|), on the order of a million years, but 
still much shorter than the relevant time-scale for plan- 
etary instability that we find §6.11 Hence from the view 
of planetary instability in situ formation and migration 
both happen close to instantaneously. 

6. INSTABILITY OF A HIGH-MASS SECOND POPULATION 

A possible way to get a moderately or highly inclined 
population of massive planets, as discussed in the previ- 
ous section, is through the intrinsic instability of plan- 
etary systems. If planetary systems are born with a 
range of characteristic masses, depending for example 
on the mass or metallicity of the protoplanetary disk, 
then the systems of higher mass can be pushed over the 
instability limit for planet-planet interaction. Planet- 
planet interactions have been studied in a wide range 
of conte xts, e.g. by intrinsic instability after gas ac- 
cretio n ([Weidenschilling fe Marzarll 119961 IRasio fc Fordl 
I1996D. in connection w ith scattering of planetesimal belts 
(|Ravmond et al.H2009D and after perturbation by a pass- 
ing s tar (jZakamska fc Tremaind 120041 : iMalmberg et al.l 
120111) . 

Long-term integration of triple systems shows that the 
separation in terms of mutual Hill radii is an impor- 
tant par ameter to determine the stability of planetar y 
systems ([Chambers et alJll996t iSmith fc Lissauerll2009| ). 
The mutual Hill radius of two planets of mass Mi and 
Mi and semi-major axis a\ and a,2 can be defined as 

\ 3M* J 2 v ' 

We choose to define the Hill radius at the mean of the 
planetary semi-major axes, following convention in the 
literature, although we note that this definition docs 
not converge towards the Hill radius of a single massive 
planet orbited by a massless test particle in the limit 
AI2 — > 0. The motion of two planets on circular or- 
bits is stable to their mutual gravity when the relative 
separation A = (a 2 - ai)/i? H > 2^3 (|Gladman|[l99l . 
Triple-planet systems do not display a similar abrupt 
transition from instability to stability. The logarithm 
of the time-scale for instability increases approximately 
linearly, log(£/Gyr) w bA + c, as the mean mutual sepa- 
ration A of the inner and outer planet p airs is increased. 
For sy stems of three low-mass planets, IChambers et al.l 
(| 19961 ) find b 1.2 and c « -1.7. Stability over at least 
10 9 years then implies A > 9. 

In Fig. [9] we show the mutual separation of planet 
neighbours, in units of their mutual Hill radii, for the 
K epler triple-transit systems (after applying our selection 
criteria). The masses are based on a simple mass- radius 
relationship proposed bv lTremaine fc Pond ()201lD to fit 
the planets in the solar system as well as the transiting 



exoplancts with known masses, 

log(i?/i?j) = 0.087 + 0.141 log(M/Mj) 

-0.171 (log M/Alj) 2 . (11) 

The radius-mass relationship increases monotonically up 
to planets of around the mass of Jupiter, and then turns 
over with a relatively constant radius for higher-mass 
planets. The high-radius branch can not be easily in- 
verted to obtain the mass, but the triple-planet systems 
in the Kepler data contain very few planets above 10 
Earth radii, so the degenerate mass-radius relationship 
is not an issue for our calculations. All systems in Fig. [9] 
are stable or at the edge of stability over 10 9 years when 
using their nominal masses. The maximum mean Hill 
separation is approximately 45, but 89% of the systems 
have (A) < 30. The lack of systems of high Hill separa- 
tion is likely partially an observational bias, since widely 
spaced systems are less likely to be observed in triple 
transit. 

Higher-mass versions of the Kepler triple-planet tem- 
plates would, with a sufficient mass boost, be unstable 
to planet-planet scattering. These systems would not 
show up in the triple-transit systems observed by the 
Kepler mission because they have lost one or more plan- 
ets by ejection or collisions (between planets or with 
the host star), leaving behind excited remnants of the 
original systems. We investigate the stability of mass- 
boosted triple-transit Kepler systems by performing TV- 
body si mulations using t he orbital dynamics code MER- 
CURY (|Chamberslll99l . We define the mass boost as 
MB = Mi/Mio where Mio is the mass of planet i from 
an approximate mass-radius relationship and Mi is the 
boosted mass. The MERCURY code uses a symplectic 
integrator to achieve conservation of linear and angular 
momentum and is thus an optimal tool for evolving plan- 
etary systems for a high number of orbital periods. 

6.1. Mass-Boosted Systems 

The triple-planet systems observed by the Kepler mis- 
sion have been carefully checked for dyn amical instabil- 
ity using the nominal plane tary masses (|Lissauer et al.l 
I2011bt iFabrvckv et"a!1l2012[ ). Here we report results of 
orbital integration of mass-boosted counterparts to the 
observed triple-planet systems. In Fig. [TO] we plot the 
results of TV-body simulations of four representative sys- 
tems from the template triple-planet systems - KOI-156, 
KOI-757, KOI-829 and KOI-408. These systems are in- 
dicated with orange dots in Fig. M We run each system 
for a number of mass boosts and monitor the time for 
the first close encounter between two planets. We con- 
sider this the relevant time-scale for instability of the 
planetary system. Each system is run for 40 random 
initialisations of the planetary orbits. The different im- 
plementations show instability at different times, but the 
median time-scale increases monotonously with the mean 
Hill separation between planet neighbours, (A). We stop 
the simulations either after 10 8 years or after 10 9 years. 
Results can be extrapolated to 10 10 years, which we take 
as an upper limit for the age of the Kepler target stars. 
Although the extrapolation has a la rge scatter between 
the in dividual systems, as noted by iDuncan fc Lissauerl 
(fl997l) who performed similar mass-boosted simulations 
of the Uranian moon system, the overall statistical trend 



Can Planetary Instability Explain the Kepler Dichotomy? 



7 



for the instability time-scale to increase rapidly with in- 
creased Hill separation is robust. 

Fig. [10] shows that the systems need (A) < 6 in or- 
der to become unstable within 10 10 years. Th is limit 
is much lower than found by IChambers et al.l (|1996f ) 
from integrations of triple systems of objects with masses 
compara ble to planetary embryo s (M w 0.03M®). As 
noted in IChambers et al.1 ()1996l ). increasing the char- 
acteristic planetary mass, while maintaining the mean 
Hill separation, leads to systematically longer instability 
time-scales. In Fig. [TT] we show the intrinsic masses of 
the planets in the triple-planet template systems, based 
on the approximate mass-radius relationship in equa- 
tion (fTTj) . as a function of the mean of the semi-major 
axis. We also indicate the masses of the systems when 
boosted to expected instability at a mean Hill separa- 
tion of (A) 6. The mass-boosted systems have char- 
acteristic planetary masses between 0.1 and 10 Jupiter 
masses, with the majority of the planets between 0.3 and 
3 Jupiter masses. 

There arc three possible (non-exclusive) outcomes 
when a planetary system becomes unstable: 1) the sys- 
tem may eject one or more planets; 2) two or more plan- 
ets may collide; or 3) one or more planets may be scat- 
tered into the host star. The relative frequency of these 
three outcomes is a function of the ratio of the planetary 
orbital speeds to their surface escape speeds, the square 
of which is often referred to as the Safronov number ( e . g . 
iSafronov fc ZviaginalH96llBinnev k Tremaind[2008l) 



Note that this defin ition is the inverse of the one used in 
IFord k R asio (2007|. In our definition, for planetary sys- 
tems having 9 > 1, planet-planet scatterings can result 
in the ejection of planets. However for systems where 
8 < 1, collisions will be common as a significant de- 
flection would only occur for close encounters (between 
two point masses) with minimum separations smaller 
than the planetary radii. The surface escape speed for 
a Jupiter-like planet is about 60 km/s, whereas the or- 
bital speed at 1 AU around a solar-mass star is 32 km/s; 
and at 0.1 AU about 100 km/s. Hence unstable plane- 
tary systems containing planets of Jovian mass will eject 
planets if their orbits have semi-major axes around or 
above 1 AU, whilst systems containing planets or tighter 
orbits around 0.1 AU will undergo planetary collisions 
(leading to mergers) or have planets scattered into their 
host stars. Collisions are thus expected to be the most 
common event for the mass-boosted planets in the Kepler 
data. 

In Fig. [T^] we show the result of integrating KOI-408 
with a mass boost of 150 (corresponding to a mean Hill 
separation (A) = 4) . The bottom panel shows the semi- 
major axes of the planets after planetary instability has 
relaxed the system to a stable configuration. The x-axis 
shows the results of forty initial representations of the 
mass-boosted system, sorted by increasing time for in- 
stability. The most common result is that the inner and 
middle planets collide. Only in systems having a long 
instability time-scale do we find collisions involving all 
three planets. The middle and top plots show the asso- 
ciated eccentricities and mutual inclinations. Eccentrici- 



ties up to 0.3 are common, while the mutual inclination 
varies from a few degrees to up to 30 degrees. This way 
the typical inclination is increased significantly by planet- 
planet scattering, so that the resulting double-planet sys- 
tems produce many observed single transits. High ec- 
centricities are also common amon g giant exoplanets ob- 
served in radial velocity surveys (jUdrv k Santos! |2007|) 
and our results show that planet-planet scattering in situ 
in the inner AU can cause both high eccentricities and 
inclinations. 

Fig.lTBIshows the fate of planets in all four dynamically 
evolved systems as a function of the mass boost (and the 
eccentricity and inclination for selected systems). Green 
denotes planet-planet collisions, blue denotes planet-star 
collisions and red ejections. The most common events 
are planet-planet collisions, as expected for these close- 
in systems. Ejections are far less common, but their fre- 
quency increases with increasing mass boost. Collisions 
with the star are very rare. 

The result of planet-planet scattering is most often 
a collision between two planets, reducing the number 
of planets from three to two. Taking this result as a 
recipe, we can create the double-planet systems arising 
from dynamical instability of triple systems by removing 
a planet from the template triple-planet systems. The 
probabilities for individual planets to transit remain un- 
changed since the (small) effect of the planetary radius 
is ignored when calculating the transit probability, al- 
though the number of single-transit, double-transit and 
triple-transit systems change when a planet is removed 
(the latter number obviously falls to zero) . The resulting 
transit numbers are shown in Fig. 1141 We experiment by 
removing either the middle planet or the outer or inner 
planet. For simplicity we maintain the original planetary 
orbits, although a more detailed analysis should put a 
planet resulting from a collision on an intermediate or- 
bit between the two colliding planets. The dashed lines 
in Fig. 03] furthermore show the result of removing the 
smallest planet in the systems. Triple transits vanish in 
all these two-planet systems. As the inclination param- 
eter P is increased the number of single transit-transit 
systems increases at the cost of double-transit systems. 

The removal of a random planet in each system would 
reduce the number of single-transit systems to 2/3 and 
the number of double-transit systems to 1/3 compared 
to the template triple-planet systems. Thus we expect 
that the ratio of single-transit to double-transit sys- 
tems will approximately double compared to the full 
systems, independently of (3. Fig. [lj] shows that it is 
possible to observe 600 single-transit systems, with only 
10 double-transit systems accompanying, if /? is moder- 
ately high at around 40 degrees. This way the massive 
double-planet systems resulting from planetary instabil- 
ity of triple-planet systems can produce a large num- 
ber of single-transit systems with massive planets and at 
the same time avoid polluting the double-transit systems 
with planets from the massive population. 

7. SYSTEMS WITH A DISTRIBUTION OF MULTIPLICITIES 

The considerations in the previous section have an in- 
herent problem in that the majority of the missing plan- 
ets in single-transit systems are small (smaller than four 
Earth radii) and thus would not contribute to making 
a system unstable. Fig. [TTJ shows that only systems 
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with characteristic planetary masses between 0.1 and 10 
Jupiter masses can undergo dynamical instability. Plan- 
etary instability of triple-planet systems of high birth 
masses can in principle result in a population of double- 
planet systems that produce the 250 observed single- 
transit systems containing large planets (more than four 
times the radius of the Earth), if the high- mass popu- 
lation represented approximately 50% of the number of 
low-mass triple-planet systems that have survived intact 
until today. Fig. [9] shows that there is a void of plane- 
tary systems with mean mutual Hill separation (A) < 9. 
This corner may have been occupied after planet forma- 
tion and become unstable to planet-planet scattering as 
the gaseous protoplanctary disk dissipated. Extrapolat- 
ing the number of systems between (A) = 9 and (A) = 18 
and between (A) = 27 and (A) = 36 to (A) < 9 shows 
that a significant number of systems could have formed 
there. 

However, the high-mass systems can not explain the 
surplus of single-transit systems with smaller planets. 
Double-transit systems and triple-transit systems also 
display a general lack of high-mass planets (see Fig. [2]), 
a property that is difficult to reconcile with the hypoth- 
esis that some systems are simply born with high-mass 
planets, since triple-planet systems could have one or two 
massive planets and still be stable over long time-scales. 

A route to producing the single-transit systems with 
small planets without making too many double-transit 
and triple-transit systems (which are already produced 
from the triple-planet system templates) is to con- 
sider additional components of double-planet and single- 
planet systems, together with the triple-planet tem- 
plates. Choosing a random inclination parameter even 
smaller than (3 = 5° for the triple-planet systems, the 
number of double-transit systems arising from the triple- 
planet templates can be reduced. Additional populations 
of double-planet and single-planet systems must then 
be present to produce the remaining single-transit and 
double-transit systems. In Fig. [15] we show the resulting 
system architectures. We make the assumption that all 
systems have the same mutual inclination parameter, a 
reasonable assumption since we have already shown that 
the double-transit and triple-transit systems in the Ke- 
pler data are statistically similar with regards to plane- 
tary sizes and orbits. We can then reconstruct the intrin- 
sic number of single, double and triple systems by requir- 
ing first that the triple-planet systems must produce all 
the triple-transit systems, then that the double-planet 
systems must produce the double-transit systems that 
are not already produced by the triple-planet systems, 
and finally that the single-planet systems must produce 
the single-transit systems not produced by double-planet 
and triple-planet systems. We create sparser systems 
from the triple-planet templates by removing one or two 
random planets from each system. 

Fig. [T5] shows that as the mutual inclination parame- 
ter is varied from 0° to 4° , the number of triple-planet 
systems increases from 3,000 to over 8,000. At the 
same time their contribution to double-transit and single- 
transit systems increases. This happens at the expense 
of double-planet systems whose number falls from ap- 
proximately 5,000 to less than zero at (3 = 5°. The 
number of single-planet systems increases from slightly 
above 20,000 to almost 27,000 as j3 increases from to 



4°. This implies that a significant fraction, more than 
20%, of the 156,000 Kepler target stars have detectable 
planets within 0.5 AU, and that the majority of planet- 
hosting stars have only a single planet in this region. We 
caution, however, that this number relies strongly on the 
assumption that rich systems of moderate or high mutual 
inclination do not contribute significantly to the observed 
transits. 

Systems of higher multiplicity than three are ignored 
in the above analysis. Without performing full synthetic 
transit observations on such systems, which is undesir- 
able since their number and statistical significance is low, 
we can derive useful limits to their contribution to triple, 
double and single transits by considering the case (5 = 0°. 
Completely flat systems have the property that the prob- 
ability that system j (of multiplicity /) displays exactly 
i transits is 

Pij = — for i = I, (13) 

P ij = ^±-^2- for i<I. (14) 

<Hj a {i+i)j 

The transit of the outermost planet always implies that 
all I planets transit (equation [T3|) . while the transit of 
planet i (with i < I) corresponds to exactly i planets 
transiting only if planet i + 1 does not transit (equation 

Analysing this way the 14 quadruple-transit systems 
which are present in the Kepler data after our selec- 
tion criteria, we find that the intrinsic population of 
quadruple-planet systems is around 400, far less than 
the 3,000 triple-planet systems present at f3 = 0° (sec 
Fig. [T5"|) . These 400 quadruple-planet systems pro- 
duce 14 quadruple transits, 8 triple-transits, 16 double 
transits and 30 single-transits, using the two transit- 
rate equations above. A similar analysis on the 4 
quintuple-transit systems yields an underlying popula- 
tion of 130 quintuple-planet systems, which give rise to 
2.3 quadruple-transit systems, 3.4 triple-transit systems, 
4.9 double-transit systems and 9.3 single transit systems. 
While systems of higher multiplicity than three must be 
present in the Kepler target stars, we can conclude that 
if the planetary systems are intrinsically very fiat, then 
their contribution to the total number of planetary sys- 
tems as well as to the lower-order transit systems is low. 
However, if the intrinsic mutual inclinations are non-zero, 
then a large fraction of the observed triple-transit sys- 
tems may arise from more crowded planetary systems 
(jRagozzine fc Holmanl[2010l ). 

Planetary systems can nevertheless not be made ar- 
bitrarily flat, unless the eccentricities in the system are 
equally low. Secular oscillations cause exchange between 
eccentricity and inclination and hence do not allow very 
flat systems. Still, the need for very flat triple-planet 
systems (as well as some systems of higher multiplicity) , 
with additional populations of low-mass single-planet 
and double-planet systems is appealing in that the abun- 
dance of planetary systems is continuous in the multiplic- 
ity. 

8. SUMMARY AND CONCLUSIONS 

This paper sets out to construct underlying planetary 
system architectures that explain the systems of single 



Can Planetary Instability Explain the Kepler Dichotomy? 



9 



and multiple transiting planets in the Kepler data. One 
of the most robust features of the Kepler data is that 
there is a dichotomy between the planets in systems dis- 
playing two or three transits and the planets in systems 
displaying only one transit. Double-transit and triple- 
transit systems have mostly small (less than four Earth 
radii) planets, while the single-transit systems host both 
small planets and an additional population of large plan- 
ets. 

The dichotomy between large planets in single systems 
and small planets in systems with a wide range of multi- 
pliciti es was also det e cted in the tra nsit synthesis mod- 
els oflLTssauer et"aLl (|2011bl) and of iTremaine fc Pond 
(|20T1 . but we confirm it here using an alternative ap- 
proach to constructing synthetic planetary systems. In 
this approach we assume that the systems of three tran- 
siting planets observed by the Kepler mission can serve 
as templates for physical planetary systems. We thus 
work close to the original data and do not have to go 
through steps of interpretation and analysis of the sta- 
tistical properties of planet sizes and orbits to fit these 
with distribution functions. 

We furthermore assume that systems with three tran- 
sits inherently arise from physical triple-planet sys- 
tems. This assumption may in reality not hold, since 
quadruple-planet systems and higher contribute to triple 
transits. However, we argue that it is reasonable to as- 
sume that denser systems are rare and thus only make 
minor contributions to observed triple-transit systems. 
The sa me assumption also exclu des the isotropic solu- 
tions of ITremaine fc Pond (|2011l) where planets are ei- 
ther in single systems or in extremely packed systems 
with almost isotropic orbits. 

We explore the underlying hypothesis that triple- 
planet systems form with a range of characteristic masses 
and that the systems of high mass are inherently unsta- 
ble. Using A-body simulations we find that the most 
likely outcome of planetary instability is collision be- 
tween two of the planets. Collisions are more likely than 
ejections because the planetary orbits arc so close to the 
host stars. The resulting double-planet systems have 
moderate mutual inclinations and would produce a num- 
ber of single transits from large planets. We find that the 
instability time-scale is very long unless we use high mass 
boosts. Triple-planet systems would need to have char- 
acteristic planetary masses between 0.1 and 10 Jupiter 
masses in order to produce dynamical instability within 
10 billion years. 

Another generic problem with invoking a second pop- 
ulation of massive planets to explain the observed large 
planets in single-transit systems is that the observed 



double-transit and triple-transit systems generally lack 
planets above four Earth radii, while these systems could 
remain stable with one or even two gas-giant planets. We 
consider this a major weakness for the hypothesis that 
the large planets in single-transit systems are a result 
of planetary instability of the high-end tail of planetary 
birth sizes. 

We find instead that the dichotomy can be seen as 
the contrast between a continuous population of almost 
flat systems with one or more small planets within 0.5 
AU versus a population of inherently single systems with 
large planets. In this picture the excess of large plan- 
ets in single-transit systems may arise already at the 
planet formation stage, when the formation or migration 
of a massive gas giant in a system suppresses the forma- 
tio n of additional (smal l and large) planets (as proposed 
by ILatham et all |20TTI) . This scenario is sketched in 
Fig. 1161 Such a mechanism has been put forward to have 
frustrated the growth of Mars in the solar system dur- 
ing J upiter's period of inwards migration (|Walsh et al.l 
l20ll . but may have been even more effective in extra- 
solar systems where gas giants have migrated to sub-AU 
orbits. 

The Kepler dichotomy appears to be ultimately con- 
trolled by the mctallicity of the host star, as gas- 
giant planets in sub-AU orbits arc mainly present 
above a threshold metallicity around the solar value 
(ISantos et alj|2004t IPTscher fc Valentll 120051: ISousa et al.l 
120 111) , while Neptune-sized planets, super-Earths and 
terrestrial planets are found in the Kep ler data around 
stars with a wide range of metallicities (|Buchhave et al.l 
120121) . 
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Fig. 1. — Triple-transit, double-transit and single-transit systems in the 16-months Kepler data after applying our selection criteria. 
Selection criteria for maximum orbital period and planet radius are indicated with dashed lines. The third selection criterion for signal-to- 
noise ratio generally excludes small planets. Systems reduced in planet number are marked with a yellow plus, while filled circles indicate 
systems that had originally the same number of planets (blue circles) or more planets (red circles). 




Fig. 2. — Planetary radius versus orbital period of the selected triple-transit, double-transit and single-transit systems. Planets in double- 
transit and triple-transit systems are statistically similar, while planets in single-transit systems are on the average significantly larger than 
in more crowded systems. 



Can Planetary Instability Explain the Kepler Dichotomy? 13 

0.251 I 




i 



Fig. 3. — The probability distribution function of the mutual inclination between pairs of planetary orbits inclined by a random angle 
distributed evenly between and /3 = 5° and with a random longitude of ascending node. The absolute maximum relative inclination 
is i = 2/3, obtained for a planet pair whose orbits are both inclined by the angle /3 but with opposite longitude of ascending node. The 
distribution peaks at a mutual inclination around i ?s 3.5° . 




Fig. 4. — The probability for the selected triple-planet systems to be seen in triple transit, for two values of the mutual inclination 
parameter /3. The system with the highest triple-transit probability, approximately 11% at /3 = 0°, is KOI-1835, a relatively large star 
orbited by planets in 2.2, 4.6 and 6.8 day orbits. The system with the lowest transit probability has a more than ten times lower probability 
for /3 = 0° and would hence be 10 times more frequent among the synthetic systems. The probability contrast increases by more than a 
factor ten when increasing the mutual inclination parameter to f) = 5°, as planets in wide orbits are much more likely to be inclined out 
of view than planets in tight orbits. 
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Fig. 5. — The number of double-transit systems from the synthetic planet population versus the number of single- transit systems, divided 
by the number of triple-transit systems, for different values of the mutual inclination parameter 0. The same ratio in the Kepler data, 
after applying selection criteria, is indicated by a star. It is not possible for any value of to get simultaneously the observed fraction of 
double-transit and single-transit systems. 
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Fig. 6. — Cumulative distribution functions of planetary radii (top) and semi-major axes (bottom) of synthetic transits (red) and of 
Kepler transits (blue) in single-transit systems (left panel), double-transit systems (middle panel) and triple-transit systems (right panel). 
We have chosen a mutual inclination parameter p = 5° that approximately matches the ratio of double to triple transits in the Kepler 
data. The synthetic transits match the radii and semi-major axes of Kepler transits by definition for triple-transit systems. Planetary radii 
show an equally good match for double- transit systems, while the single-transit systems have significantly larger planets in the Kepler data. 
Single-transit and double-transit systems have observed planets slightly further from the star than the synthetic planets. 
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Fig. 7. — Cumulative distribution function (top panel) and probability distribution function (bottom panel) of the absolute number 
of observed planets (blue) and synthetic planets (red) in single-transit systems versus their radius. The difference between observed 
and synthetic transits is shown in yellow. Approximately 750 single-transit planets are missing in the synthetic population, which only 
reproduces 1/3 of the actual number of planets in single-transit systems. Approximately 500 of the missing planets are small (smaller than 
four Earth radii) while 250 are large (larger than four Earth radii). 
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Fig. 8. — Number of synthetic single-transit, double-transit and triple-transit systems versus mutual inclination parameter /3, normalised 
to match the number of observed triple-transit systems at f} = 5°. The isotropic limits are indicated with dashed lines. Triple-transit 
systems are quickly suppressed when increasing /3, but almost isotropic inclinations are needed to suppress double-transit systems. Hence 
is is difficult for a second population of highly inclined triple-planet systems to avoid overproducing (massive) double-transit systems. 
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Fig. 9. — The separation, measured in mutual Hill radii, of planet neighbours in triple-planet systems. The triple-planet stability criterion 
for Earth-mass planets - a mean Hill separation of 9 — is given by a dashed red line for the generic planet mass based on an approximate 
mass-radius relationship. Mutual Hill radii of [18,27,36,45], corresponding to (A) = 9 for mass boosts of [8,27,64,81], are also indicated. 
The systems that are evolved dynamically in J6]are shown in orange. 
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Fig. 10. — Time-scale for a close encounter between two planets to occur in mass-boosted triple-planet systems, as a function of the mean 
initial separation in units of mutual Hill radii, of four representative triple-planet systems (KOI-156, KOI-757, KOI-829 and KOI-408). 
Dots show the results of different initial realisations of the planet orbits for a given (A), with big crosses indicating the median instability 
time-scale and pluses the mean (the bar shows the standard deviation). Individual mass boosts are printed next to the mean, while the 
top axis shows the mass boost for KOI-156. The number of implementations of each system that are stable for more than either 10 8 or 
10 years are shown at the top of the plot next to an upwards-pointing arrow. Most systems need a mean Hill radius separation (A) < 6 
to become unstable within 10 10 years. 
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Fig. 11. — Masses of planets in the triple-planet template systems as functions of the mean semi-major axis of planets in the system. 
The blue points show the intrinsic masses based on an approximate mass-radius relationship, while the red points show the masses after 
boosting each system to expected instability at a mean Hill separation of (A) = 6. The boosted planetary masses are between 0.1 and 10 
Jupiter masses, with the majority of the planets between 0.3 and 3 Jupiter masses. 
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Fig. 12. — Result of orbital integration of KOI-408. The bottom panel shows the semi-major axis after planetary instability (left axis, 
planets indicated with circle radius proportional to planetary radius) and the time-scale for the first ejection or collision (right axis, turquoise 
curve). The x-axis shows the results of forty initial representations of the mass-boosted system, sorted by increasing time for instability. 
An empty column at position 31 is used to separate systems with two remaining planets from systems with one remaining planet. Typically 
the middle planet collides with the inner planet, or in cases where the system falls apart slowly, all three planets merge. The middle and 
top panels show the eccentricities and mutual inclinations of the post-encounter systems. The resulting eccentricity is typically ~0.1, while 
the mutual inclination ranges between 1 and 30 degrees. 
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FlG. 13.— Outcome of planetary instability of KOI-156, KOI-757, KOI-408 and KOI-829 for different mass boosts. KOI-156 was 
furthermore run for four different values of the eccentricity and three values for the inclination and KOI-829 for two values of the eccentricity. 
Each parameter choice was run for 40 different random initialisations. Green denotes planet-planet collisions, blue denotes planet-star 
collisions and red/brown ejections. Planet- planet collisions are by far the most common event due to the proximity to the star. Ejections 
are less common, but their frequency increase with increasing mass boost. Collisions with the star are very rare. 
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Fig. 14. — The transit numbers after removing a planet by planet-planet collision, collision with the star or ejection. The blue lines 
show the number of single-transit systems, when the middle (full line) or inner/outer planet (thin line) is removed. The red lines show the 
number of double-transit systems. The dashed line shows the result of removing the smallest planet in the system. Triple-transit systems 
vanish. The number of double-transit systems is reduced, compared to using the original triple-planet systems, relatively more than the 
number of single-transit systems. 
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Fig. 15. — Fractional contribution to transits (left panel) and number of planetary systems (right panel), as functions of mutual inclination 
parameter (horizontal axis) and system multiplicity (yellow bars: single-planet systems, red bars: double-planet systems, blue bars: 
triple-planet systems). Here we have relaxed the assumption that all planetary systems are triple and allowed for both single-planet, 
double-planet, and triple-planet systems. Increasing 0, the number of triple-planet systems increases, as does their contribution to single 
and double transits. The number of double-planet systems decreases slowly as increases, eventually dropping below zero for = 5° (not 
shown). The necessary population of large single-planet systems is marked with a dotted line. 
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Fig. 16. — Sketch of the various situations that can lead to a separate population of massive single-planet systems. In the top panel triple- 
planet systems form with various masses, either in situ (left) or further out followed by migration (right). The more massive systems are 
unstable and reduce the planet number by planet-planet collisions, leaving two planets of moderately high mutual inclination. However, the 
masses needed to make the systems unstable are very high. Instead the formation or migration of a gas-giant planet may suppress formation 
of other (small and large) planets in the sub-AU regions of the system, since the gas giant excites high eccentricities and inclinations of the 
planetcsimal and embryo population (lower panel). 



